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Theory  and  Simulation  of  the  Rayleigh-Taylor 
Instability  in  the  Large  Larmor  Radius  Limit 

I .  INTRODUCTION 

Interchange  instabilities,  such  as  the  Rayleigh-Taylor  instability, 
have  been  the  subject  of  intense  scientific  interest  for  most  of  this 
century. 1  Generically,  an  interchange  instability  can  occur  when  two 
fluids  of  different  densities  are  accelerated.  Although  initially  studied 
in  the  context  of  hydrodynamic  theory, 2  this  type  of  instability  has  also 
been  shown  to  exist  in  magnetized  plasmas.-^  Subsequently,  plasma 
interchange  modes  have  become  an  important  area  of  research  in  the  plasma 
physics  community,  and  a  considerable  amount  of  work  (theoretical, 
experimental,  and  computational)  has  been  done  over  the  past  two  decades. 
The  bulk  of  the  research  has  been  motivated  by  the  fusion  research  program, 
especially  in  regard  to  the  stability  of  magnetically  confined  plasmas^*  and 
imploding  laser  target  pellets.-*  However,  there  is  also  evidence  that  this 
type  of  instability  can  occur  in  space  plasmas.^  Thus,  interchange  modes 
are  relevant  to  a  broad  spectrum  of  plasma  regimes. 

The  study  of  the  Rayleigh-Taylor  instability  in  a  magnetized  plasma  has 
largely  been  confined  to  the  regime  pj/L  <<  1  and  <<  1  where  and 

are  the  ion  gyroradius  and  gyrof requency ,  respectively,  and  w  and  L  are  the 
frequency  and  length  scales  of  interest,  respectively.  Standard  MHD  theory 
can  be  used  when  kp^  <<  1  (where  k  is  the  wavenumber)  while  f ini te-Larmor- 
radius  (FLR)  MHD  theory  is  used  when  kpj  <  1.  [Of  course,  use  of  kinetic 
theory  can  remove  all  restrictions  on  the  value  of  kp^.®’^]  However, 
space*®  and  laboratory** ’ *^  experiments  have  observed  plasma  structure, 
apparently  caused  by  the  rapid  deceleration  of  a  plasma  expanding  into  a 
strong  magnetic  field,  in  the  parameter  regime  Pj/L  >>  1  and  w/ 52^  >>  1  but 
Pe/L  <<  1  and  «/S2e  <<  1.  Similar  structure  has  also  been  observed  in 

earlier  plasma  expansion  experiments*-^’*^  and  in  theta  pinch  implosion 
experiments*-*  in  the  regime  p^/L  <  1  and  w/S^  <  1.  The  purpose  of  this 

Letter  is  to  investigate,  both  analytically  and  computationally,  the 
development  of  the  Rayleigh-Taylor  instability  in  the  large  ion  Larmor 
radius  regime  (i.e.,  p^/L  >>  1  and  w/S?^  >>  1).  We  show  that  in  this  regime 
the  Rayleigh-Taylor  instability  develops  much  faster  and  has  a  dramatically 
different  nonlinear  behavior  than  the  conventional  Rayleigh-Taylor 
instability. 
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II.  THEORY 

A  set  of  one-fluid,  modified  MHD  equations  has  been  derived  which 
describes  a  plasma  for  arbitrary  p^/L  or  w/Qj : ^ 

|y  +  V  •  nu  =  0  (1) 

du  j 

nM  dT  =  ~  T7n  “  73  /811  +  8  '  7B/4n  +  nMg  (2) 

3B  M  du 

aF  =  9  x  x  fi)  -  r  v  x  af  (3) 

where  u  is  the  center  of  mass  velocity,  n  is  the  density,  B  is  the  magnetic 
field,  M  is  the  ion  mass,  and  g  is  the  gravitational  force.  [With  regard 
to  plasma  expansion  experiments,  we  can  identify  the  plasma  acceleration 
with  an  effective  gravitational  acceleration,  i.e.,  du0/dt  =  -  geffl- 
Equations  (1)  -  (3)  are  derived  directly  from  the  ion  and  electron  Vlasov 
equations  under  the  following  assumptions.  We  assume  (1)  pe/L  <<  1  and 
«/Qe  <<  1  (i.e.,  the  electrons  are  strongly  magnetized),  (2)  the  electron 
and  ion  pressure  tensors  are  isotropic,  (3)  the  electron  and  ion  fluids  are 

isothermal  with  T  =  Te  +  Tj,  and  (4)  the  Debye  length  is  very  small.  We 

mention  that  (1)  -  (3)  are  similar  to  a  set  of  one-fluid  MHD  equations 
derived  in  Ref.  18  in  the  limit  Tj  =  0. 

The  final  term  in  (3)  is  what  makes  these  equations  different  from  the 

conventional  ideal  MHD  equations.  In  this  representation  it  is  clear  that 
ion  inertia  is  the  origin  of  this  term  and  that  it  is  potentially  important 
when  w/Qj  >>  1  (or  p^/L  >>  1).  We  also  point  out  that  (3)  can  be  rewritten 
as 


8B 

;rr  =  V  x  (u  x  B)  -  7  x  (|  x  B)/ne 


so  that  the  new  term  can  be  identified  as  the  Hall  term. 


(4) 
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To  study  the  linear  evolution  of  the  Rayleigh-Taylor  instability  in 
both  the  small  and  large  Larrnor  radius  limits,  we  consider  the  following 
slab  geometry  and  plasma  configuration.  Ve  assume  B  =  Bq(x)  ez,  n  =  ng(x), 
and  g  =  g  ex.  The  equilibrium  is  given  by  3/3x  (ngT  +  Bq^/8h)  =  ngMg.  Ve 
perturb  (1)  -  (3)  about  this  equilibrium  and  assume  perturbations  are 
proportional  to  exp  [i(kyy  -  cot)].  We  consider  the  "local"  limit  and 
require  kyLn  >>  1  where  Ln  =  (31n  ng/3x)“^  is  the  scale  length  of  the 
density  gradient.  The  dispersion  equation  is  given  by^6 


* (i  *  «“2  -  tyt1  - 

*  r  rr  (“2  -  ky«Ln)“  ■ 0 


where  =  B/(4itnM)^^^  is  the  Alfven  velocity,  Cs  =  (2T/M)^^  is  the  sound 
speed,  pj  =  Cs/Sj,  and  3  =  8llnT/B^. 

In  the  limit  pj/Ln  -»  0  (i.e.,  conventional  MHD  theory)  we  recover  from 

(5)  the  magnetosonic  modes  [co^  =  kyV^(l  +  3)]  and  the  usual  Rayleigh-Taylor 
instability  [co^  =  g/Ln(l  -  31nBg/ 31nng) ] .  Instability  can  occur  when  g/Ln 
<  0  with  the  growth  rate  yo  =  (g/Ln)^^.  However,  when  pj/Ln  »  1,  a 

faster  growing  mode  can  occur  which  has  a  growth  rate  determined  by  co^  = 
kygLn.ll> Again,  instability  can  occur  when  gLn  <  0  but  the  growth  rate 
is  y  =  ky(gLn)l^  =  kyLpYgi  the  turn-on  condition  for  this  mode  is 
(g/Ln)  >  Interestingly,  for  a  magnetically  confined 

plasma  expansion  in  a  low  3  plasma,  the  turn-on  condition  is  independent  of 
the  magnetic  field  which  is  consistent  with  experimental  observations.^ 
Of  further  interest  is  the  eigenmode  structure.  We  note  that  when  pj/Ln  >> 
1  and  co  ~  kyCs,  then  w/Sj  >>  1  and  the  final  term  in  (3)  will  be  large 

unless  ez  •  V  x  u  ->  0.  We  use  this  fact  in  the  equation  obtained  by  taking 
the  z  component  of  the  curl  of  (2)  to  get  -  ioo(  3n0/3x)Uy  =  ikygn.  From  (1) 
we  obtain  -  icon  +  ingkyUy  +  (3ng/3x)ux  =  0.  But  since  z  •  V  x  u  =  0  we 

note  that  ikyUx  =  0  so  that  -  icon  +  ingkyUy  =  0.  Thus,  eliminating  n  and 
Uy  from  these  equations  we  get  co^  =  kygLn.  The  crucial  point  is  that  the 
fluid  motion  is  primarily  curl-free  and  not  divergence- free  (i.e.,  V  x  u  = 
0  and  7  •  u  =  ikyUy  *  0).  This  represents  a  behavior  quite  different  from 

the  conventional  MHD  Rayleigh-Taylor  instability  which  is  primarily 
divergence-free  and  not  curl-free  (V  •  u  =  0  and  V  x  u  ?  0). 


This  contrasting  fluid  behavior  can  also  be  seen  by  considering  the 
following,  simplified  ion  momentum  equation  in  the  limit  T  -»  0:  nMdVj/dt  = 
e(E  +  Yi  x  B/c).  In  the  latter  case,  when  w/Qj  <<  1,  the  dominant 
perturbed  ion  velocity  is  =  -cE  x  B/c  which  is  divergence-free.  In  the 
former  case,  when  w/SJj  >>  1,  the  dominant  perturbed  ion  velocity  is  = 
ieE/nMw  which  is  curl-free.  Thus,  based  upon  this  linear  eigenfunction  we 
anticipate  that  the  nonlinear  behavior  of  the  unmagnetized  and  magnetized 
ion  Rayleigh-Taylor  instabilities  will  be  quite  different. 

We  now  confirm  this  suspicion  by  presenting  results  of  2D  numerical 
simulations  of  the  Rayleigh-Taylor  instability  in  both  the  small  Larmor 
radius  limit  (i.e.,  magnetized  ions)  and  the  large  Larmor  radius  limit 
(i.e.,  unmagnetized  ions).  We  solve  (1),  (2),  and  (4)  numerically  on  a 
cartesian  grid  with  resolution  125  x  100  (x,y).  We  consider  plasma  motion 
only  in  the  plane  transverse  to  the  magnetic  field  (B  =  -Bez).  Hard  wall 
boundary  conditions  were  used  in  the  x  direction,  while  periodic  boundary 
conditions  were  used  in  the  y-direction.  The  algorithm  used  solves  (1), 
(2),  and  (4)  in  conservation  form  using  a  Total  Variation  Decreasing  (TVD) 
non-linear  switch  between  a  first  order  transport  scheme  and  an  Adams- 
Bashforth,  centered  8th  order  spatial  scheme. ^  The  following  density 
profile  was  chosen  ng(x)  =  0.55  +  0.45tanh(x/L2)  for  x  <  0,  and  nQ(x)  = 
[0.55  +  0.45tanh(x/L2) lexpC-x/Lj)  for  x  >  0,  with  g  =  -  g  ex  and  B(x) 
determined  from  equilibrium  force  balance.  We  choose  =  1.0,  L2  =  0.2, 
Cs  =  1.0,  and  g  =  1.0.  The  length  scale  is  normalized  to  Lj  =  1.0,  the 
velocity  is  normalized  to  Cs  =  1.0,  and  the  time  scale  is  normalized  to 
Li/Cs  =  1.0.  The  boundaries  used  in  the  simulations  are  x  =  -2.5  and  2.5, 
and  y  =  0.0  and  4.0. 

We  benchmarked  the  code  in  the  linear  regime  by  initializing  the  grid 
with  a  single  mode  and  calculating  the  eigenfunction  and  eigenvalue.  For 
Pi/L2  =  0  and  Ax  =  .8  we  calculate  the  linear  growth  rate  to  be  y  ~  1.70 
which  agrees  well  with  the  theoretical  value  of  y  ~  1.81.  In  the  nonlinear 
regime  it  is  difficult  to  benchmark  the  code  with  previous  results  because 
the  code  is  compressible  and  we  treat  a  diffuse  density  profile  (as  opposed 
to  a  sharp  boundary  model).  Nevertheless,  we  estimate  the  Froude  number  to 
be  F  ~  0.27  which  is  in  reasonable  agreement  with  experimental^'  (F  ~  0.2  - 
0.3)  and  numerical  (F  ~  0.23)  results.  The  initial  conditions  for  the 
simulations  discussed  in  the  remainder  of  the  Letter  is  a  1%  random  density 
fluctuation  throughout  the  grid. 
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In  Fig.  1  we  show  the  density  contours  (Fig.  la)  and  velocity  vector 
field  (Fig.  lb)  for  the  Rayleigh-Taylor  instability  in  the  small  Larmor 
radius  limit  (i.e.,  we  use  Pj/L^  =  0)  at  t  =  5.5.  We  show  these  results 
primarily  for  comparison  with  the  Rayleigh-Taylor  instability  in  the  large 
Larmor  radius  limit.  The  following  points  should  be  noted.  First,  the 
instability  initially  develops  in  the  region  x  <  0.2.  The  peak  of  the 
density  at  t  =  0  is  at  x  =  0.2  so  that  the  linearly  unstable  regime  is  x  < 
0.2  where  g31nnQ/3x  <  0.  As  the  instability  evolves  nonlinearly,  we  find 
that  the  usual  'spike'  and  'bubble'  morphology  develops  (see  Fig.  la).^2 
The  'heavy'  material  falls  (i.e.,  the  'spike')  while  the  'light'  material 
rises  (i.e.,  the  'bubble').  In  each  case,  a  mushroom-like  structure 
develops.  Second,  associated  with  the  development  of  the  mushroom-like 
heads  is  the  development  of  counter-rotating  vortices.  This  is  evident  in 
Fig.  2a.  The  fluid  motion  in  the  'spikes'  is  such  that  the  fluid  falls 
downward  in  the  center  of  the  spikes,  but  at  the  edges  is  deflected 
horizontally,  and  eventually  turns  upwards  creating  a  swirling  pattern. 
The  'bubbles'  display  the  same  behavior  but  the  upward  flowing  fluid  is 
eventually  turned  around  to  go  downward.  Thus,  the  fluid  has  a  significant 
curl  as  expected.  And  third,  the  bulk  of  the  'action'  takes  place  in  the 
linearly  unstable  region  (i.e.,  x  <  0.2).  We  see  that  disturbed  fluid 
motion  occurs  between  x  =  -  1.0  and  x  =  0.7  so  by  this  time  there  is 
relatively  little  disruption  of  the  'stable'  side  of  the  density  profile 
(i.e.,  where  g31nnQ/3x  >  0). 

In  Fig.  2  we  show  the  density  contours  (Fig.  2a)  and  velocity  vector 
field  (Fig.  2b)  for  the  Rayleigh-Taylor  instability  in  the  large  Larmor 
radius  limit  (i.e.,  we  use  p^/L2  =  5.0)  at  t  =  1.0.  The  following  points 
should  be  noted.  First,  it  is  clear  that  the  unmagnetized  Rayleigh-Taylor 
instability  develops  much  faster  than  the  magnetized  Rayleigh-Taylor 
instability,  consistent  with  linear  theory.  Second,  although  not  evident 
in  Fig.  2,  we  find  that  during  the  linear  evolution  of  the  unmagnetized 
instability,  the  fastest  growing  modes  are  the  shortest  wavelength  modes. 
This  is  consistent  with  the  theoretical  analysis  which  indicated  that  y  a 
ky.  However,  longer  wavelength  structures  seem  to  dominate  as  the 
instability  progresses,  i.e.,  those  that  are  comparable  to  the  initial 
density  gradient  scale  length.  Third,  the  instability  initially  develops 
on  the  'unstable'  side  of  the  density  profile  (i.e.,  g31nnQ/3x  <  0),  but 
unlike  the  magnetized  instability,  it  disrupts  more  of  the  density  profile 


(see  Fig.  2a).  The  extent  of  the  disruption  is  between  x  =  -  1.1  to  x  = 
1.6,  a  significantly  larger  region  than  the  previous  case.  It  appears  that 
narrow  channels  develop  which  allow  the  magnetic  field  to  penetrate  into 
the  'stable'  topside  of  the  density  profile.  Fourth,  a  secondary 
instability  appears  to  develop  on  one  side  of  the  largest  density 
enhancements,  i.e.,  the  'fishbone-like'  structure  in  Fig.  2a.  It  is  not 
entirely  clear  what  cause  this  finer  scale  structure;  it  may  be  an 
eigenmode  structure  or  a  secondary  instability.  Fifth,  the  velocity  vector 
field  is  completely  different  from  the  magnetized  ion  case.  The  flow  is 
clearly  not  divergence-free  (see  Fig.  2b)  and  large-scale  flows  seem  to 
form  fluid  'sheets'.  The  maximum  velocity  reached  in  this  case  is  um  =  1.9 
which  is  somewhat  larger  than  the  'free-fall'  velocity  u  =  1.5.  And 
finally,  we  add  that  we  have  also  performed  simulations  for  p}/L2  ~  0(1); 
the  linear  growth  rate  of  the  instability  is  smaller  than  either  of  two 
examples  presented  above  in  accordance  with  linear  theory. ^ 

III.  CONCLUSION 

In  conclusion,  we  have  presented  the  theory  and  first  nonlinear  fluid 
simulations  of  the  Ray leigh-Taylor  instability  in  the  large  ion  Larmor 
radius  limit  (i.e.,  p1/L  >>  1).  These  results  are  based  upon  a  modified 
set  of  one-fluid  MHD  equations  which  is  valid  for  arbitrary  values  of  pj/L 
or  <0/9}  (but  requires  pe/L  <<  1  and  <o/9e  <<  1).  We  have  demonstrated  that 
in  this  limit  the  Ray leigh-Taylor  instability  evolves  in  a  dramatically 
different  manner  from  the  conventional,  small  Larmor  radius  Rayleigh-Tay lor 
instability.  Besides  developing  much  faster,  the  unmagnetized  instability 
appears  to  be  more  disruptive  of  an  accelerating  plasma  density  shell,  and 
generates  compressible  ion  flows.  Moreover,  although  the  fastest  growing 
modes  have  the  shortest  wavelength,  longer  wavelength  modes  seem  to 
dominate  in  the  nonlinear  regime.  Finally,  these  results  are  consistent 
with  the  observed  structuring  of  the  AMPTE  barium  release  in  the 
magnetotail , ^  and  recent  laboratory  observations  of  plasma  structure.^ 
In  both  cases,  plasma  shells  rapidly  decelerate  as  they  expand  into  a 
strong  magnetic  field,  develop  structure  on  time  scales  much  faster  than 
the  ion  cylcotron  period,  and  are  dominated  by  structure  sizes  comparable 
to  the  density  gradient  scale  length. 
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Fig.  2)  Simulation  results  for  the  Ray leigh-Taylor  instability  in  the 
large  Larmor  radius  limit  (pj/L^  =  5)  at  time  t  =1.0.  The 
magnetic  field  is  in  the  -7.  direction  (out  of  the  page).  (a) 
Density  contours.  The  contour  levels  are  the  same  as  in  Fig.  la. 
(b)  Velocity  vector  field.  The  'unit'  velocity  vector  of  1.5  is 
shown  for  comparison. 
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0 1 C  Y  ATTN  DONALD  HANSEN 


INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA  ,  VA  22  31  1 
0 1 C  Y  ATTN  J.M.  AEIN 
01 CY  ATTN  ERNEST  BAUER 
0 1 C  Y  ATTN  HANS  WOLFARD 
0 1 C  Y  ATTN  JOEL  BENCSTON 

I  NT  L  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY ,  NJ  071 1 0 

0 1 C  Y  ATTN  TECHNICAL  LIBRARY 

JA  YCOR 

11011  TORREYANA  ROAD 

P.O.  BOX  85159 

SAN  DIEGO,  CA  921 38 

0 1 C  Y  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
0 1 C  Y  ATTN  THOMAS  POTEMRA 
0 1 C  Y  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7963 

COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

8 1 6  STATE  STREET  (P.O  DRAWER  QQ) 
SANTA  BARBARA,  CA  93102 
01CY  ATTN  DAS  I A  C 
01CY  ATTN  WARREN  S.  KNAPP 
01CY  ATTN  WILLIAM  MCNAMARA 
01CY  ATI N  B.  CAMBILL 

LINKABIT  CORP 
10953  ROSELLE 
SAN  DIEGO,  CA  92121 

01CY  ATTN  IRWIN  JACOBS 
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LOCKHEED  MISSILES  &  SPACE  CO.,  INC 

P.0.  BOX  50*1 

SUNN  XV  Al.F,  ,  CA  94088 

0  1  C  Y  ATTN  DEPT  60-12 
0  1  C  Y  ATTN  D.H.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3  2  5  1  HANOVER  STREET 
PALO  ALTO,  CA  9430^ 

01CY  ATTN  MARTIN  WALT  DEPT  52-12 
0 1 C  Y  ATTN  W.L.  IMHOF  DEPT  52-12 
01 CY  ATTN  RICHARD  G.  JOHNSON 
DEPT  52-12 

01CY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  6 1  37 
ORLANDO,  FL  32805 

0 1 C  Y  ATTN  R.  HEFFNER 

MCDONNEI.  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  9264  7 
01CY  ATTN  N.  HARRIS 
0  1  C  Y  ATTN  ,1.  MO II I, E 
OILY  ATTN  GEORGE  MROZ 
01CY  ATTN  W.  OLSON 
0 1 C  Y  ATTN  R.W.  HALPRIN 
01 CY  ATTN  TECHNICAL 

LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
0 1 C  Y  ATTN  P.  FISCHER 
0 1 C  Y  ATTN  W.E.  CREVIER 
0 1  C Y  ATTN  STEVEN  !,.  GUT.SCHE 
01CY  ATTN  R.  BOGUSCH 
0 1 C  Y  ATTN  R.  HENDRICK 
01CY  ATTN  RALPH  KILB 
0 1 C  Y  ATTN  DAVE  SOWLE 
01  CY  ATTN  F.  FA  JEN 
0 1 C  Y  ATTN  M.  SCHEIBE 
0 1 C  Y  ATTN  CONRAD  L.  LONGMIRE 
0 1 C  Y  ATTN  B.  WHITE 
01 CY  ATTN  H.  ST A  GAT 

MISSION  RESEARCH  CORP. 

1  /'MO  K  A  N  D  0 ! ,  P  H  ROAD,  S  .  E  . 

ALBUQUERQUE,  N  M  H  7  /  0 o 
(.'ICY  h.  M.TELL  I  NCWE  RE 

oil'  y  .  a;, me 

0  1  C  Y  WRIGHT 


MITRE  CORP 

WESTCATE  RESEARCH  PARK 
1 820  DOLLY  MADISON  BLVD 
MCLEAN ,  VA  22101 

01CY  ATTN  W.  HALL 
01 CY  ATTN  W.  FOSTER 

PACIFIC -SIERRA  RESEARCH  CORP 
1  2  3*10  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01 CY  ATTN  K.C.  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  l.AB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVEHSiTY  PAHK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01 CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRIES,  INC. 
a  ARROW  DRIVE 
WOBURN,  MA  01001 

01  C  Y  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

0 1 C  Y  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC  . 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

0 1 C  Y  ATTN  FORREST  GILMORE 

0  1 C  Y  ATTN  WILLIAM  B.  WRIGHT,  JR. 

01CY  ATTN  WILLIAM  J.  KARZAS 

0 1 C  Y  ATTN  H.  ORY 

0 1 C  Y  ATTN  C.  MACDONALD 

0 1 C  Y  ATTN  BRIAN  LAMB 

0 1 C  Y  ATTN  MORCAN  GROVER 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY  ,  MA  01776 

01CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  IN.  ilTUTE 
3  ?  0  W  E  S  T  4  2  n  :i  2,  r  f<  E  E  T 
NEW  Y  5RK  ,  NY  1 0036 

0  1  C  Y  ATTN  VI  M C E  T R A  P A  N  I 


SCIENCE  APPLICATIONS 

INTERNATIONAL  INCORPORATED 
1 1 50  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

0 1  C Y  ATTN  LEWIS  M.  LINSON 
0 1 C  Y  ATTN  DANIEL  A.  HAMLIN 
0 1 C  Y  ATTN  E.  FRIEMAN 
0 1 C  Y  ATTN  E . A .  STRAKER 
0 1 C  Y  ATTN  CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS 

'INTERNATIONAL  CORPORATION 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
0 1 C  Y  J.  COCKAYNE 
0 1 C  Y  E.  HYMAN 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK  ,  CA  99  025 


01  C  Y 

ATTN 

J  . 

CASPER 

01CY 

ATTN 

DONALD  NEILSON 

01  CY 

ATTN 

ALAN  BURNS 

01CY 

ATTN 

G  . 

SMITH 

01  CY 

ATTN 

R. 

TSUNOLA 

0 1  C  Y 

ATTN 

DAVID  A.  JOHNSON 

01  CY 

ATTN 

WALTER  G.  CHESNUT 

0 1  C  Y 

ATTN 

CHARLES  L.  RINO 

0 1  C  Y 

ATTN 

WALTER  JAYE 

0  1  C  Y 

ATTN 

J  . 

VICKREY 

01  CY 

ATTN 

RAY 

L.  LEADABRAND 

0 1  C  Y 

ATTN 

G  . 

CARPENTER 

01  CY 

ATTN 

G  . 

PRICE 

0  1  C  Y 

ATTN 

R. 

LIVINGSTON 

01  CY 

ATTN 

V  . 

GONZALES 

0  1  C  Y 

ATTN 

D  . 

MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD  ,  MA  01730 

0 1 C  Y  ATTN  W.P.  BOQUIST 

TRW  DEFENSE  A  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
n  E  DO.\  DO  BEACH,  CA  90278 


0  1  C  Y 

ATTN 

R  . 

K.  PLEBUCH 

01  CY 

ATTN 

S  . 

ALTSCHULER 

0  1  C  Y 

ATTN 

D  . 

DEE 

01  CY 

ATTN 

D/ 

STOCKWF.I.L 

SNTP /I  5  7  5 

V  IS  I  DYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MA  01803 
01 C  Y  ATTN  W .  RE  I DY 
0 1 C  Y  ATTN  J.  CARPENTER 
01 C  Y  ATTN  C  .  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH ,  PA  15213 

0  1  C  Y  ATTN  :  N  .  ?.  ABUSKY 


